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Abstract — We propose a kinematic approach to searching for the stars that could 
be formed with the Sun in a common "parent" open cluster. The approach consists 
in preselecting suitable candidates by the closeness of their space velocities to the solar 
velocity and analyzing the parameters of their encounters with the solar orbit in the past 
in a time interval comparable to the lifetime of stars. We consider stars from the Hipparcos 
catalog with available radial velocities. The Galactic orbits of stars have been constructed 
in the Allen-Santillan potential by taking into account the perturbations from the spiral 
density wave. We show that two stars, HIP 87382 and HIP 47399, are of considerable 
interest in our problem. Their orbits oscillate near the solar orbit with an amplitude of 
~250 pc; there are short-term close encounters to distances < 10 pc. Both stars have an 
evolutionary status and metallicity similar to the solar ones. 



INTRODUCTION 

According to present-day observations, stars are formed in groups or clusters of various 
strengths. It is highly likely that the Sun was formed in some "parent" open star cluster. 

Comparison of the solar chemical composition with the metallicity distribution in the 
Galactic disk led to the conclusion that the Sun was born ^2 kpc closer to the Galactic 
center relative to its current position (Wielen et al. 1996). A discussion of this problem 
in light of the currently available data can be found in Acharova et al. (2010). Revealing 
possible siblings of the Sun will allow a better understanding of both the mechanisms of 
radial migration in the Galactic disk (Minchev and Famaey 2010; Shevchenko 2010) and 
the conditions under which the Sun was formed (Williams 2010). It is hypothesized that 
the stars born together retain their chemical homogeneity for a very long time (Bland- 
Hawthorn and Freeman 2004; Bland-Hawthorn et al. 2010). 

Therefore, as a method of searching for such stars, these authors suggested analyzing 
spectroscopic data to determine the abundances of various elements. In particular, the 



spectroscopic method has been successfully used to prove the dynamical origin of the stars 
belonging to the Hyades cluster (Pompeia ia et al. 2011). 

Obviously, the stars formed together with the Sun are also of great interest in studying 
both the dynamical evolution of the Solar system and the development of terrestrial life 
(Valtonen et al. 2009). 

The difficulties of revealing such stars in the solar neighborhood are associated with 
the Sun's great age (p^4.6 Gyr), the dynamical evolution of open star clusters (OSCs) 
during this time, and a large uncertainty in Galactic parameters. Numerical simulations 
of the dynamical evolution of OSCs show (Chumak et al. 2005; Chumak and Rastorguev 
2006a, 2006b) that stellar tails stretched along the Galactic cluster orbits develop in them 
with time. The OSC remnants existing in the form of tails must completely dissolve and 
mix with the stellar background in a time ~2 Gyr (Kiipper et al. 2008). According to 
the estimates by Portegies Zwart (2009), 10-60 stars from the parent (for the Sun) OSC 
containing ps 10 3 members can now be in a solar neighborhood about 100 pc in radius. 

The attempts at searching for the Sun's siblings from a common cluster are known in 
the literature. For example, Brown et al. (2010) searched for suitable candidates among 
the stars of the Hipparcos catalog (1997) in an updated version (van Leeuwen 2007). A 
list of candidates containing six stars was proposed. However, when the orbital motion 
of stars was simulated, the influence of the perturbations from the spiral density wave 
was disregarded. In addition, the search was carried out using only the stellar parallaxes 
and proper motions (without radial velocities). Mishurov and Acharova (2010) showed 
that the influence of the spiral density wave in a time interval of 4.6 Gyr leads to a 
significant dispersal of the members of an initially compact cluster in both radial and 
tangential directions. For about a hundred dispersed stars to be now observable in a solar 
neighborhood p^IOO pc in radius, the parent cluster must contain ps 10 4 members. 

In this paper, we suggest applying a kinematic approach to searching for the stars 
that could be formed with the Sun in a common parent cluster. This approach consists 
in (1) selecting suitable candidates by the closeness of their space velocities to the solar 
velocity, (2) constructing the stellar and solar orbits in the Galactic potential including 
the perturbations from the spiral density wave, and (3) analyzing the parameters of the 
encounters between the stellar orbits and the solar orbit in the past in a time interval 
comparable to the lifetime of stars. 

DATA 

The initial set of kinematic data on the stars from the Hipparcos catalog (1997) that 
we used is described in detail in Bobylev et al. (2010). It contains the parallaxes and 
proper motions of about 35000 stars taken from a revised version of the Hipparcos catalog 
(van Leeuwen 2007). The radial velocities from the PCRV catalog (Gontcharov 2006) are 
available for each star. 

We selected 162 F, G, and K stars with a relative parallax error <j n /n < 15% and mag- 
nitude of the total stellar space velocity (U, V, W) relative to the Sun y/U 2 + V 2 + W 2 < 
8 km s^ 1 , where 8 km s _1 was estimated from a typical random error in each of the 
velocities U, V, W, which are p^2 km s" 1 . Note also that age and metallicity estimates are 
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available for a significant fraction of the selected stars (Holmberg et al. 2009). 

Remarkably, two of the six candidates from Brown et al. (2010) entered into our list 
of 162 stars: HIP 21158 and HIP 99689. 



THE CONSTRUCTION OF ORBITS 

We calculated the stellar and solar orbits by solving the following system of equations of 
motion based on a realistic model of the Galactic gravitational potential (Fernandez et 
al. 2008): 

£ = --^-^(i?o -0-2^077, (1) 
i) = + + 2% 

where $ is the Galactic gravitational potential; the (£, i], () coordinate system with the 
center at the Sun rotates around the Galactic center with a constant angular velocity 
Q , with the £,77, and ( axes being directed toward the Galactic center, in the direc- 
tion of Galactic rotation, and toward the Galactic North Pole, respectively; R is the 
Galactocentric distance of the Sun. 

We used the Allen-Santillan (1991) model Galactic potential. The system of equa- 
tions (1) was solved numerically by the fourth-order Runge-Kutta method. 

In the Allen-Santillan (1991) model, the Galactocentric distance of the Sun is taken 
to be R = 8.5 kpc and the circular velocity of the Sun around the Galactic center is 
Vq = \Qq\Rq = 220 km s -1 . The axisymmetric Galactic potential is represented as the 
sum of three components — the central (bulge), disk, and halo ones: 

$ = $ c + $d + $ff- (2) 

The central component of the Galactic potential in cylindrical coordinates (r, 9, z) is 
represented as 

$ c = M£ (3) 



where M c is the mass and be is the scale parameter. The disk component is 
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where is the mass, and &d are the scale parameters. The halo component is 

M g (it>/q g ) 2 - 02 
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Here, is the mass and an is the scale parameter. If R is measured in kpc and 
Mq, Md, Mh are measured in units of the Galactic mass (Mq) equal to 2.32 x 10 7 M Q , 
then the Gravitational constant is G — 1 and the unit of measurement of the potential 
$, along with the individual components of (3)-(5), is 100 km 2 s~ 2 . 

All of the Allen-Santillan model parameters adopted here are given in Table 1. 

If the spiral density wave is taken into account (Lin and Shu 1964; Lin et al. 1969), 
then the following term is added to the right-hand side of Eq. (2) (Fernandez et al. 2008): 

$ sp (R,6,t) = Acos[m(n p t - 6) + x(R)}, (6) 

where 

_ ( J R ^o) 2 / r0 tani 

A — ) 
m 

, ^ m / R\ 

Here, A is the amplitude of the spiral wave potential; f r Q is the ratio of the radial compo- 
nent of the perturbation from the spiral arms to the Galaxy's total attraction; Q p is the 
pattern speed of the wave; m is the number of spiral arms; % is the arm pitch angle, i < 
for a winding pattern; x is the phase of the radial wave (the arm center then corresponds 
to x = 0°); and x& is the Sun's phase in the spiral wave. 

The spiral wave parameters are very unreliable (a review of the problem can be found 
in Fernandez et al. (2001) and Gerhard (2010)). The simplest model of a two-armed 
spiral pattern is commonly used, although, as analysis of the spatial distribution of young 
Galactic objects (young stars, star-forming regions, or hydrogen clouds) shows, both three- 
and four-armed patterns are possible (Russeil 2003; Englmaier et al. 2008; Hou et al. 
2009). More complex models are also known, for example, the kinematic model by Lepine 
et al. (2001) that combines two- and four-armed spiral patterns rigidly rotating with 
an angular velocity close f2 - Note also the spiral ring Galactic model (Mel'nik and 
Rautiainen 2009) that includes two outer rings elongated perpendicular and parallel to 
the bar, an inner ring elongated parallel to the bar, and two small fragments of spiral 
arms. As applied to the Galaxy, the theories of nonstationary spiral waves with a fairly 
short stationarity time (several 100 Myr), a variable rotation rate, and a variable number 
of arms are also considered (Sellwood and Binney 2002; Baba et al. 2009). The currently 



available data do not yet allow one of the listed models to be unequivocally chosen. 
Therefore, here we apply the model of a stationary spiral pattern with different numbers 
of spiral arms. 

Given the data on the Galactic bar rotation (Debattista et al. 2002), the pattern speed 
Hp can lie within the range 15-65 km s _1 kpc -1 . The possibility of the coexistence of 
several pattern speeds, a rapidly rotating bar and a slower spiral pattern, is also considered 
(Minchev and Famaey 2010; Gerhard 2010). Here, we disregard the influence of the bar. 
Therefore, we choose fl p from the range 15-30 km s^ 1 kpc -1 (Popova and Loktin 2005; 
Naoz and Shaviv 2007; Gerhard 2010). 

The pitch angle % is known relatively well and is —5° -= — 7° and —10° -= — 14° for the 
two- and four-armed spiral patterns, respectively. 

The amplitudes of the velocities of the perturbation from the spiral density wave are 
5 — 10 km s" 1 (Mishurov and Zenina 1999; Fernandez et al. 2001; Bobylev and Bajkova 
2010) in both tangential and radial directions. 

The Sun's phase in the wave x& is known with a very large uncertainty. For example, 
according to Fernandez et al. (2001), this angle lies within the range 284° — 380°. Having 
analyzed the kinematics of OSCs (Bobylev et al. 2008; Bobylev and Bajkova 2010), 
we found x& close to —117° (or 243°), while the data on masers yielded an estimate of 
Xq = —130 ± 10°. (Bobylev and Bajkova 2010); here, the minus implies that hat we 
measure the phase angle from the center of the Carina-Sagittarius spiral arm. From the 
position of the Sun among the spiral arms, it is obvious that this angle is close to — n/2 
(Russeil 2003). 

According to the classical approach in the linear density- wave theory (Yuan 1969), the 
ratio f r o lies within the range 0.04-0.07 and the most probable value is f r o = 0.05. The 
upper limit f r o = 0.07 is determined by the velocity dispersion of young objects observed 
in the Galaxy (at f r o = 0.07, the dispersion must reach 25 km s _1 , which exceeds a typical 
observed value of 10-15 km s" 1 ). 

We adopted the following parameters: the two-armed spiral pattern (m = 2), the 
pattern speed of the spiral wave Q p = 20 km s _1 kpc -1 , the pitch angle % = —5°, the 
ratio f r0 = 0.05, and the Sun's phase in the wave x® — —117°. We also took into 
account the Sun's displacement from the Galactic plane Z Q = 17 pc (Joshi 2007) and 
used the present-day peculiar velocity of the Sun relative to the local standard of rest 
([/©, V Q , W q ) lsr = (10,11,7) km s^ 1 (Binney 2010; Schdnrich et al. 2010; Bobylev and 
Bajkova 2010). 

For all of the selected 162 stars, we determined the relative distances d between the 
star and the Sun as well as their velocity difference dV in the time interval in the past 
with a boundary of —4.5 Gyr. The minimum values of these parameters, d min and dV min , 
were also determined for a certain time t min . 

RESULTS 

Encounter Parameters 

Figure 1 presents the (B — V) color-absolute magnitude My diagram for the 162 selected 
candidates. To calculate My, we used the apparent V magnitudes from the Hipparcos 
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Note. For the age, the upper and lower limits of the estimate calculated with a 1 error are 
given in parentheses; 1 — Holmberg et al. (2009), where the metallicity was determined from 
Stromgren photometry; 2 — Robinson et al. (2007), where the metallicity was determined from 
Lick indices; 3— Takeda et al. (2007). 



catalog. We see from Fig. 1 that the number of candidates lying close to the isochrones 
and having nearly solar spectral types will be considerably smaller than the total number 
of stars selected by their kinematics at the first step. 

The distributions of the sample of 162 stars in the Galactic XY, ZX, and YZ coor- 
dinate planes are shown in Fig. 2. The X, Y, and Z coordinate axes are directed toward 
the Galactic center, in the direction of Galactic rotation, and toward the Galactic Pole, 
respectively. 

We found that only two stars from our sample, namely HIP 47399 (Fig. 3) and 
HIP 87382 (Fig. 4), can be acceptable candidates. For them, there is good agreement 
between their estimated ages and the kinematic encounter parameters (d, dV) found. In 
Figs. 1 and 2, the positions of these stars are marked by the large filled circles (the random 
error bars are within the circles) and the crosses, respectively. 

For HIP 47399, we have the following initial data: a = 9 /l 39 m 27 s .4, 5 = 42°17'09", 
fi a cos5 = —4.29 ± 1.38 mas yr" 1 , fis = —6.70 ± 0.59 mas yr _1 , n = 13.87 ± 0.95 mas 
(the heliocentric distance r = 72 ± 5 pc), the radial velocity V r = —7.3 ± 3.4 km s _1 , 
U = 3.9 ± 2.3 km s _1 , V = -2.5 ± 0.3 km s _1 , W = -6.3 ± 2.6 km s _1 . 

For HIP 87382: a = 17 fe 51 m 14 s .O, 5 = 40°04'20", /i a cos5 = -16.86 ± 0.28 mas yr" 1 , 
Us = 11.01 ± 0.42 mas yr" 1 , vr = 29.76 ± 0.36 mas (r = 33.6 ± 0.4 pc), Vr = 1.6 ± 0.2 km 
s -1 , U = -1.0 ± 0.1 km s -1 , V = 0.4 ± 0.2 km s~\ W = 3.4 ± 0.1 km s -1 . The ages and 
metallicities of these stars estimated by various authors are given in Table 2. The error 
in [Fe/H] is typically 0.1-0.15 dex. This suggests that the metallicity of HIP 47399 and 
HIP 87382 is nearly solar, within the error limits. 

HIP 21158 and HIP 99689 common to the list by Brown et al. (2010) do not withstand 
the described test for close encounters in the past. 

The encounter parameters d and dV depend on the adopted components of the Suns 
peculiar velocity vector ([/©, V & , Wq)lsr an d on whether or not we take into account the 
influence of the spiral wave. Without allowance for the spiral wave, there are no close 
encounters between HIP 47399 and the Sun, while the closest encounters for HIP 87382 
{dmin = 16 pc, dVmin = 16 km s _1 ) occur at tmin = —1-2 Gyr. It may be concluded that 
allowance for the spiral density wave has a significant influence on the results. 
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Figure 1: Positions of stars on the colorabsolute magnitude diagram; the zero- age main 
sequence (4) and three isochrones (Demarque et al. 2004) for ages of 3 (3), 4 (2), and 5 
(1) Gyr are marked. 



The encounter parameters calculated for HIP 47399 and HIP 87382 by taking into 
account the influence of the spiral wave are given in the last columns of Table 2. 

The Probability of Close Encounters 

To estimate the probability that a star was born together with the Sun in a common 
star cluster, we should take into account the following. The lifetime of the OSC as a 
gravitationally bound structure is typically no more than 2 Gyr. The OSC core size is 
typically about 10 pc. The OSC has an open halo and a tail elongated along its orbit 
develops with time. We assumed that the OSC was born at a time of —4.5 Gyr. 

We then consider the probability (p) that the star and the Sun belong to the OSC 
core to be nonzero if the following conditions are met for them in the time interval 
—4.5 < t < —2.5 Gyr: the relative distance along the Galactic radius is \dR\ < 20 pc, 
along the vertical axis is \dZ\ < 20 pc, and along the orbit is \dY\ < 20 pc, then d = 
V dR 2 + dZ 2 + dY 2 < 35 pc. For HIP 87382, according to the upper panel in Fig. 4, 
these conditions are close to fulfilment in this time interval — we observe a minimum of 
d = 40 — 60 pc near —3.5 Gyr. Similar reasoning is also valid for HIP 47399. 

As a result, we calculate p as the ratio of the number of points in the narrow peaks 
(Figs. 3 and 4) satisfying the adopted constraints to the total number of points in a given 
time interval. 

The constraint on the relative velocity at the encounter time should also be taken 
into account. Then, for example, for HIP 87382 at d < 10 pc and \dV\ < 10 km s^ 1 
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Figure 2: Distributions of stars in the XY, ZX, and planes; the filled square and 
cross mark the stars HIP 47399 and HIP 87382, respectively. 
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Figure 3: Parameters d and dV of the encounter between HIP 47399 and the solar orbit 
versus time for the two-armed spiral pattern. 



in the time interval from —4.5 to —2.5 Gyr, the probability p = 0.001(0.1%) is then a 
more rigorous estimate. For less rigorous requirements, d < 20 pc and \dV\ < 10 km s" 1 , 
p = 0.006(0.6%). 

For HIP 47399, these numbers are slightly larger (because its relative velocity is always 
low, dV < 10 km s -1 , as we see from the lower panel in Fig. 3): p — 0.006(0.6%) for 
d < 10 pc and \dV\ < 10 km s" 1 and p = 0.014(1.4%) for d < 20 pc and \dV\ < 10 km 
s" 1 . 

For a rigorous estimation of p, we should take into account the errors in the obser- 
vational data and the errors in the parameters of the Suns peculiar velocity and the 
parameters of the spiral density wave using the method of statistical simulations. In 
Fig. 5, the parameters of the mutual encounters between HIP 87382 and HIP 87382 are 
plotted against time. We can see that their relative velocity in the time interval from 
—3.5 to —2.8 Gyr is 5-10 km s -1 ; there is a peak encounter to distances of less than 10 
pc. Hence it may be concluded that within the framework of the approach used, three 
objects (the Sun, HIP 47399, and HIP 87382) are possible candidates for being members 
of a common open cluster. 
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Figure 4: Parameters d and dV of the encounter between HIP 87382 and the solar orbit 
versus time for the two-armed spiral pattern. 



Statistical Simulations 

We calculated the encounter parameters by taking into account the random errors. 
The errors in the Suns peculiar velocity components (U q ,V & ,W q )lsr were taken to be 
(0.5, 1, 0.3) km s -1 , a 10% level in the initial velocities (U, V, W) and coordinates (X, Y, Z) 
of the stars as well as in the parameters of the spiral density wave. 

We performed our simulations for HIP 87382, because it has very small random errors 
in the initial coordinates and velocities. It turned out that there are interesting results in 
our simulations even in a symmetric potential (mainly due to the variations in the Suns 
peculiar velocity). Therefore, we present the results for HIP 87382 in the time interval 
from —4.5 to —2.5 Gyr for two cases: (1) the orbits were constructed only in a symmetric 
potential and (2) with the influence of the spiral density wave added: 

1. For d < 10 pc, \dV\ < 10 km s~\ p = 0.0002(0.02%), for d < 10 pc, \dV\ < 5 km 
s" 1 , p = 0.0001(0.01%). 

2. For d < 10 pc, \dV\ < 10 km s _1 , p = 0.000009(0.0009%), for d < 10 pc, 
\dV\ < 5 km s-\ p = 0.000006(0.0006%). 

We see that p decreased by less than one order of magnitude in case 1 and by two 
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Figure 5: Parameters d and dV of the mutual encounters between HIP 47399 and 
HIP 87382 versus time for the two-armed spiral pattern. 



orders of magnitude in case 2 compared to the analogous results described in the preceding 
section for this star. Obviously, the uncertainty in the parameters of the spiral density 
wave has a decisive influence. 

One of the most important parameters in our model of the spiral wave is the ratio 
of the radial component of the perturbation from the spiral arms to the Galaxys total 
attraction, f r0 . To investigate the robustness (stability to small perturbations) of our 
results, we determined the encounter parameters for f r o from the range 0.04 — 0.07 (with 
all of the remaining model parameters being fixed). The results are presented in Table 3. 
We see from this table that the results remain acceptable for our problem up to f r0 = 0.065 
(the perturbation is 30% with respect to 0.05). 

Other Spiral-Structure Models 

The four-armed spiral pattern. For the four-armed spiral pattern, we set m = 4 
and the pitch angle % = —10°. The remaining parameters were the same as those for the 
two-armed model, namely £l p = 20 km s _1 kpc -1 , f r0 = 0.05, and x& — —117°. The 



Table 3: Encounter parameters calculated using the two-armed (m = 2) spiral pattern at 
various f r0 
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Note. * — for the local minimum from the interval t < — 2 Gyr. 

simulation results are presented in Figs. 6 and 7. 

We can seen from comparison of Figs. 6 and 3 that the four-armed model of the spiral 
pattern for HIP 47399 yields more interesting results: at t ~ 4 Gyr, a global minimum is 
observed in encounters to distances of less than 30 pc with relative encounter velocities 
< 2 km s _1 . 

As can be seen from comparison of Figs. 7 and 4, the picture of encounters with the 
solar orbit did not change qualitatively for HIP 87382. 

The 2+4 spiral pattern. For the composite (2 + 4) model of the spiral pattern 
(Lepine et al. 2001; Mishurov and Acharova 2010), the spiral wave potential (6) contains 
two terms with amplitudes A 2 (two-armed component) and A4 (four-armed component). 
In this model, the Sun is very close to the corotation circle and, as Mishurov and Acharova 
(2010) showed, the influence of the spiral structure is so strong that the test model particles 
are scattered over a significant spatial volume. 

The following model parameters were adopted. For the two-armed component: m = 
2,i 2 = —7°,x& = 300°. For the four-armed component: m = 4, ? 4 = — 14°, x = 135°. 
Both spiral patterns rotate with the same angular velocity Q p = Q . In our case, 
tt = 220/8.5 = 25.9 km s" 1 kpc^ 1 , as follows from the Allen-Santillan (1991) model 
parameters. We took f r0 = 0.05 when calculating A 2 and used A2/A4 = 0.8 to determine 
A 4 . 

The results are presented in Figs. 8 and 9. We can see that for both stars, HIP 47399 
and HIP 87382, there are encounters of interest to us. 

CONCLUSIONS 

Based on the proposed kinematic approach to searching for the Sun's siblings from a 
common cluster, we selected 162 F, G, and K stars from the Hipparcos catalog with low 
heliocentric space velocities (< 8 km s" 1 ) from a solar neighborhood ^200 pc in radius. 
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Figure 6: Parameters d and dV of the encounter between HIP 47399 and the solar orbit 
versus time for the four-armed spiral pattern. 



For all these stars, we constructed the Galactic orbits for 4.5 Gyr into the past using the 
Allen-Santillan (1991) axisymmetric model Galactic potential that additionally included 
the perturbations from the spiral density wave. The parameters of the encounter with 
the solar orbit were calculated for each orbit. 

We confirmed the conclusion reached by Mishurov and Acharova (2010) that the spiral 
density waves have a significant influence on the search results. 

We considered the influence of a stationary spiral structure for the two- and four- 
armed models m = 2,4 as well as the composite 2 + 4 model (Lepine et al. 2001) and 
found that they all confirm the results of our search. 

We found that almost all of the stars considered in the past receded from the solar 
orbit fairly rapidly and far. Two single stars, HIP 87382 and HIP 47399, constitute an 
exception. In the two-armed model of the Galactic spiral pattern, their orbits oscillate 
near the solar orbit with an amplitude of ^250 pc; there are peak close encounters to 
distances d < 10 pc; both stars have a wide minimum up to distances d < 60 pc in the 
time interval — 4-h— 3 Myr consistent with their age estimates. HIP 47399 is interesting in 
that its heliocentric velocity is always low, while this velocity in the time interval — 4-h— 3 
Myr is \dV\ < 3 km s -1 . The closest and longest encounters at t ~ —4 Gyr are observed 
for the four-armed model (m = 4) for HIP 47399 (Fig. 6). 

As a result, we conclude that HIP 47399 and HIP 87382 are of considerable interest 
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Figure 7: Parameters d and dV of the encounter between HIP 87382 and the solar orbit 
versus time for the four-armed spiral pattern. 



as possible candidates for the Sun's siblings from the hypothetical parent cluster. 
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